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HYCOM  =  HYbrid  Coordinate  Ocean  Model 
NCODA  =  Navy  Coupled  Ocean  Data  Assimilation 
=  Polar  Ice  Prediction  System 

Delivery: 

•  Scheduled  for  transition  to  the  Naval  Oceanographic  Office  in  2008 

Capability: 

•  Provide  accurate  3D  temperature,  salinity  and  current  structure 

•  Depict  the  location  of  mesoscale  features:  fronts  and  eddies 

Progress: 

•  1/1 2 “global  HYCOM/NCODA  running  in  real-time  in  the  NAVOCEANO 
operational  queues  since  22  Dec  2006 

•  Produces  daily  5-day  hindcast  up  to  the  nowcast  time,  then  a  5-7  day 
forecast 

•  Graphical  and  digital  output  available  through  the  HYCOM  consortium 
web  pages:  http://www.hvcom.org 

•  Validation  efforts  underway  -  comparing  against  operational  1/8° 
global  NCOM  (Navy  Coastal  Ocean  Model) 


Global  HYCOM  Configuration 

■  Horizontal  grid:  1/1 2 “equatorial  resolution 

4500  x  3298  grid  points,  ~6.5  km  spacing  on  average,  ~3.5  km  at  pole 

■  Mercator  79  °S  to  47  °N,  then  Arctic  dipole  patch 

■  Vertical  coordinate  surfaces:  32  for  cr2* 

■  KPP  mixed  layer  model 

■  Thermodynamic  (energy  loan)  sea-ice  model  -  switching  to  PIPS 

■  Surface  forcing:  NMOC  NOGAPS  0.5°wind  stress,  wind  speed, 
thermal  forcing,  and  NOGAPS  1.0 “precipitation 

■  Monthly  river  runoff  (986  rivers) 

■  Initialized  from  January  climatology  (GDEM3)  T  and  S,  then  SSS 
relaxation  from  PHC  3.0 

No  subsurface  relaxation  to  climatology 


Validation  Tasks 

A.  Large  scale  circulation  features 

•  Determine  correct  placement  of  large  scale  features 

B.  Sea  Surface  Height  (SSH)  variability  /  Eddy  Kinetic 
Energy  (EKE) 

•  Determine  if  the  system  has  a  realistic  level  and  distribution  of  energy 
at  depths 

C.  Mixed  layer  depth  (MLD  /  sonic  layer  depth  (SLD)  / 
deep  sound  channel  (DSC)  /  below  layer  gradient  (BLG) 

•  Compare  simulated  vs.  observed  for  non-assimilated  buoys 

D.  Vertical  profiles  of  T&S 

•  Quantitative  comparison  of  simulated  vs.  observed  for  non-assimilated 
buoys 

E.  Sea  surface  temperature 

•  Evaluate  whether  the  models  are  producing  acceptable  nowcasts  and 
forecasts  of  sea  surface  temperature 

F.  Coastal  sea  level 

•  Assess  the  model’s  ability  to  represent  observed  sea  surface  heights 


Large  Scale  Circulation  Features 


1992-2002  Mean  Dynamic  Ocean  Topography  (MDOT) 


from  Maximenko  and  Niiler  (2005) 


Large  Scale  Circulation  Features 


2004-2006  mean  sea  level  from  1/1 2°  global  HYCOM/NCODA 
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From  the  1/1 2°  global  HYCOM/NCODA  hindcast  simulation 
standard  deviation  of  difference  (MDOT  -  HYCOM)  =  9.2  cm 
standard  deviation  of  difference  (MDOT  -  NCOM)  =  13.0  cm 


SSH  Variability  Evaluation 


Measure  of  the  mesoscale  eddy  field 


Oct  92  -  May  07  SSH 
variability  based  on  T/P, 
ERS-1  and  ERS-2 
altimeters  (from 
Collecte,  Localisation, 
Satellites  (CLS)) 


SSH  variability  over 
2004-2006  from  a  1/12° 
global  HYCOM/NCODA 
hindcast  simulation 


SSH  Variability  Evaluation 


Measure  of  the  mesoscale  eddy  field 


Area  avg.  =  7.8  cm 


Area  avg.  =  6.8  cm 


Oct  92  -  May  07  SSH 
variability  based  on  T/P, 
ERS-1  and  ERS-2 
altimeters  (from 
Collecte,  Localisation, 
Satellites  (CLS)) 


SSH  variability  over 
2004-2006  from  the  1/8° 
global  NCOM  real-time 
simulation 


Surface  EKE  in  the  Gulf  Stream 


Observations  from 
Fratantoni  (2001)  - 
based  on  1 990-99 
surface  drifters 
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MLD  Error  Analysis 

JJ/2007  MLD  Median  Bias  Error  (MdBE):  model  vs.  ~6K  assim.  profiles 


HYCOM/NCODA  MdBE:  -5  m,  RMSE:  53  m 


60°N 

pO°N 

0° 

30°S 

60°S 

90°S 

60°N 

pO°N 

0° 

30°S 

60°S 

90°S 


#  of  Profiles 

- 1 - 


#  of  Profiles 


- 1 - 

60°E 


SC 


gNCOM  02f  MdBE:  -4  m,  RMSE:  62  m 
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Temp  vs.  Depth  Error  Analysis 


vs.  NCODA  analysis  in  z-space  (green) 

vs.  NCODA  analysis  in  HYCOM  space  (blue) 

Mean  error  (°C)  RMSE  (^C) 


Based  on  ~5300  assimilated  profiles  over  the  period  June-July  2007 


Sea  Surface  Temperature  Evaluation 

Data  type:  MCSST,  ~1 9,000,000  observations 


Data  type:  Drifting  buoys,  -520,000  observations 


Based  on  thirty  4-day  /  3-day  forecasts  from  HYCOM  /  NCOM  over  the  period 

June-July  2007;  Limited  between  45°S  -  45°N 


Sea  Ice  Simulation  in  HYCOM 

Existing  public  domain  version  of  HYCOM  is  configured 
with  a  thermodynamic  energy-loan  seajce  model  built  in 

■  First  generation  system 

■  No  ice  rheology  Hice  grows/melts  as  a  function  oMieat  flux  &  SST 

Couple  HYCOM/NCODA  with  a  sea  ice  (ClCij)  model 
developed  by  Los  Alamos 

■  Next  generation,  advanced  system 

■  Additional  Ice  physics 

■  Energy-based  ice  fidging  scheme 

■  Energy-conserving  thermodynamics 

■  Multi-category,  fflearly  remapped  ice  thickness 

■  2-way  coupling  between  ocean  and  ice  via  tie  Earth  System 
Modeling  Framework  (ESMF) 

■  In  Navy  parlance:  Polar  Ice  Prediction  System  (PIPS) 

■  Assimilate  SSMI  ice  concentration  in  PIPS 


Navy  Coupled  Ocean  Data  Assimilation  (NCODA) 


Ocean  Obs 


SST:  GAC/LAC 
MCSST,  GOES, 
Ship,  Buoy 
Profile:  XBT, 
CTD,  PALACE 
Float,  Fixed  Buoy, 
Drifting  Buoy 
Altimeter  SSHA 
SSM/I  Sea  Ice 


Sequential  Incremental  Update  Cycle 
Analysis-Forecast-Analysis 


Innovations 


increments 


Forecast  Fields  First  Guess 
Prediction  Errors  s 


MVOI  -  simultaneous  analysis  5  ocean  variables  temperature, 
salinity,  geopotential,  layer  pressure,  velocity  (u,v) 


HYCOM/NCODA  Runstream 

Hindcast _ _ _ Forecast 
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NCODA  analysis  windows  centered  at  these  times 
±36  hours  for  altimeter  data 
±12  hours  for  all  other  data 

1)  Perform  first  NCODA  analysis  centered  on  tau  =  -126 

2)  Run  HYCOM  for  24  hours  using  incremental  updating  (□)  over  the  first  6  hrs 

3)  Repeat  steps  1 )  and  2)  until  the  nowcast  time 

4)  Run  HYCOM  in  forecast  mode  out  to  tau  =  120 

Approximate  run  times*  (using  379  IBM  Power  5+  processors): 

1 )  Six  NCODA  analyses:  1 .1  hrs/analysis  =  6.6  hrs 

2)  Five  HYCOM  hindcast  days  @  240  sec  At:  0.8  hrs/day  =  4.0  hrs 

3)  Five  HYCOM  forecast  days  @  240  sec  At:  0.8  hrs/day  =  4.0  hrs 

4)  Total:  14.6  hrs 


*  Timings  do  not  include  PIPS  coupling 


MLD/SLD/BLG/DSC  Evaluation 


Mixed  Layer  Depth  (MLD)^ 
change  in  temperature  of 
.25qClrom  tie  surface 

Sonic  Layer  Deptf|  (SLD): 
near  surface  sound  speed 
maximum 

Below  Layer  Gradient  (BLG): 
fit  a  line  to  sound  speed 
points  between  SLD  and 
SLD  + 100  m 

Deep  Sound  Channel  (DSC) 
axis:  deepes|sound  speed 
minimum 


